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1.0 INTRODUCTION 


INPUTS is a program for interpolation in both space (2-dimensional or 3-dimensionaI) 
and time. It is based on a linear Interpolation scheme using simplex spatial regions 
(triangles and tetrahedra). INPUTS was developed mainly to provide data interfacing 
between the output from thermal analyzers and the input to the BO PACE 3-D (Boeing 
plastic analysis capability for 3-dimensional solids using isoparametric finite elements) 
program. The INPUTS interpolator is of a quite general nature, however, and could 
be used effectively for other tasks (with appropriate changes in input/output formats, 
if necessary). 

The INPUTS program utilizes temperature values which are given at some sequence of 
times for a list of strategically located "thermal nodes." It operates on these values 
by performing a double interpolation (in time and space) to provide temperature values 
at another desired sequence of times for a list of "structural nodes." The thermal/ 
structural interface procedure is shown in Figure 1-1 . 

INPUTS is written in FORTRAN IV and is available on both the IBM 360/370 and the 
UNIVAC 1108 machines. The program has a core storage requirement of 30K words, and 
it presently has a capability for handling 1000 structural nodes and 500 thermal nodes. 

The INPUTS document consists of three major parts: 

PART I Theoretical Manual 

PART II User Manual 

PART III Programmer Manual 

Recognition Is due to Curt Whiting for his work in preparation of example problems, 

and in debugging of the program. 
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Figure 1-1. Thermal/Structural Interface for INPUTB/BOPACE Run 
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2.0 INPUTB INTERPOLATION THEORY 

INPUTB performs a double interpolation in time and space> to convert data associated 
with one set of times and sfxatial locations (i.e , , thermal node data) to data associated 
with another set (i.e., structural node data). In order to progressively develop the 
necessary concepts, the timewise interpolation is discussed first, followed by^ discussions 
of the 2-dimensionat and 3-dimensIona! spatial interpolation, 

2.1 TIMEWISE INTERPOLATION 

The interpolation in time is a simple linear, 1-dimensional scheme. Given the tempera- 
ture of a point at a sequence of "thermal" times, two of these times are used to compute 
the temperature at each desired "structural" time. The two thermal times ,-of course, 
must be those nearest to the structural time, such that the first is less than or equal to, 
and the second is greater than or equal to, the structural time. 

2.2 • SPATIAL 2-DIMENSIONAL INTERPOLATION 

The 2-dimensional Interpolation scheme is illustrated by Figure 2-1 . There the open 
circles represent thermal nodes with given temperature values, and. the closed circle 
(point N) represents one of the structural nodes at which the temperature must be computed. 

In general, three thermal node points must be selected (forming an enclosing triangle 
around the structural node) in order to accomplish an exact linear interpolation. The 
sequence for selecting these three points is as follows: 

1, Find point 1 aS'the thermal node closest to N. 

2. Locate point 2 as the next closest thermal node such that the smaller angle 
between lines N-1 and N-2 is greater than 90°, (The restriction on the 
angle is necessary to avoid possible selection later of the third point such 
that a long narrow triangle would be formed. Note it also ensures that a 
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perpendicular from N to the line 1-2 will pass between points 1 and 2). 

Locate point 3 as the next closest thermal node such that N is enclosed 
within the triangle 1-2-3. (A necessary and sufficient condition for satisfying 
this requirement is that the three vector cross products N-1 X N-2, N-2 X N-3, 
and N-3 X N-1, all have the same direction) o 

The structural node N divides the triangle into three triangular sub-areas A-j, P<2, and 
A3, with total area of the triangle defined by A = A^ -^2 ^ Interpolation weighting 
factors are then defined by W] = Ai/A, W2 = A2/A, and W3 = A3/A, so that Wi + 

W2 + W3 — 1. Finally, the temperature Tjv^ at point N is computed by 

Tn = W] Tt +W 2 T2 + W3 T3 (2-1) 

where Tj is the temperature of the i th thermal node forming the triangle. 

In case a third thermal node is not found to form an enclosing triangle, a perpendicular 
is defined from N to the line 1-2, intersecting 1-2 at say point N^. The weights W-j 
and W2 are then computed as the distances from to points 2 and 1, respectively, 
divided by the length of line 1-2. Tj^ can then again be computed from Equation 2-1 
(INPUTS logic in this case sets W3 = 0, and assumes node 3 = node 1). In case a 
second node is also not found according to the above selection procedure, INPUTS 
logic sets W2 — W3 — 0, and assumes node 2 = node 3 = node 1, and again makes use of 
■ Equation 2-1 . 
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Figure 2-1. 2-Dimensional Spatial Interpolation 
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Figure 2-2. 3-Dimensional Spatial Interpolation 
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2.3 


SPATIAL 3-DIMENSlONAL INTERPOLATION 


The 3-dImensional Interpolation scheme is similar to the 2-dImensional scheme, except 
that tetrahedra rather than triangles are involved. As shown by Figure 2-2, the open 
circles again represent thermal nodes with given temperature values, and the closed 
circle (point N) represents a structural node at which the temperature must be computed. 

Four thermal node points must now be selected to form an enclosing tetrahedron around 
the structural node in order to accomplish an exact linear interpolation. The sequence for 
selecting these four points is as follows: 

1 . Find point 1 as the thermal node closest to N. 

2; Locate point 2 as the next closest thermal node such that the smaller angle between 
lines N-1 and N-2 Is greater than 90°. 

3. Locate point 3 as the next closest thermal node such that a perpendicular from N 
to the plane 1-2-3 intersects this plane within the triangle 1-2-3. (A vector S is 
defined norma! to the plane 1-2-3. The requirement is satisfied if the three vector 
dot-cross products ,S • (N-1 x N-2), S • (N-2 x N-3), and S • (N.-3 x N-1), 

all have the same sign.) 

4. Locate point 4 as the next closest thermal node such that N is enclosed within the 
tetrahedron 1-2-3. (A necessary and sufficient condition for satisfying this require- 
ment is that the four vector dot-cross products -N-1* (N-3 x N-4), + N-2*(N-4 x N-1), 
-N-3* (N-1 X N-2), and + N-4* (N-2 x N-3), all have the same sign.) 
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The structural node N divides the tetrahedron into four tetrahedral sub-volunies, V-|, 

^2' '^ 3 ' volume of the tetrahedron defined by V = V-; + V2 + Vg + 

(Vj is the volume of the tetrahedron opposite node poi nt i). interpolation weighting factor: 
are then defined by W; = V]/V, W2 = V2/V, W3 = and W4 = V4A'/ so that 

W-| + W2 + W3 + W4 = 1 . Finally the temperature at point N is computed by 

Tn = W^Ti + W2T2 + W3T3 +W4T4 (2-2) 

where Tj is the temperature of the i th thermal node forming the tetrahedron. 

In case a fourth thermal node is not found to form an enclosing tetrahedron^, point N is 
projected perpendicularly onto the 1-2-3 plane, and the 2-dimensionaI interpolation 
scheme using triangle 1-2-3 is employed. In case a third or second node is also not found, 
INPUTB logic again proceeds as discussed for 2-dimensionaI interpolation. 
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3,0 PROGRAM FLOW 


The major steps occomplished during an INPUTB run are shown tn the main program flow 
summary of Figure 3-1. Step 1 reads input/output file unit number, and allows for printed 
and/or punched output. READM and READMT read the structural and thermal mesh 
(node locations) respectively, with structural nodes given in BOPACE 3-D formot. 

Step 4 reads a series of thermal node temperature vectors, corresponding to a given 
sequence of times. In step 5, the major interpolation logic is employed, as discussed in 
Sections 2,2 and 2,3, to select four thermal nodes and compute associated weighting 
factors, for each structural node. The four nodes and weighting factors are then stored, 
so that they con later be used to give Interpolated structural node temperatures at each 
desired time. In step 6, the interpolation is actually performed, using these stored values 
and the thermal node time-temperature vectors read in Step 4, Multiple problems may be 
executed during a single program run. 
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■Another yes 

PROBLEM 

^ET? 


Figure 3-1. Main Problem Flow 
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4.0 SUMMARY OF ItMPUTB INPUT DATA 

A picforial of the INPUTS input deck is shown in Figure 4-1. The input data consists 
of the following two general types. 

Type Cl Data on the usual card file. 

Type [• Data on File 1. These include almost all' of the input data. File I 

is defined by the user so that his data can be Input via cards, tape, etc. 

The data included on these files are described below. Formats are consistent with 
FORTRAN IV conventions. 

C-1. Start code and file unit numbers: 


a. 

Insert the code "START" - 


b. 

Unit number for file 1. (Must be given). 


c. 

Unit number for first output file (e.g., printer). 

(If not given, this 


file will not be output.) 


d. 

Unit number for second output file (e.g. punch). 

(If not given, this 


file will not be output.) 


Format 

(A4, 6X, 315) 



1-1. Problem l.D. title (any 80 characters). 

Format (20A4). 

1-2 Structural node location data (BOPACE 3-D format). For each structural node 
give node l.D. number, identification number of coordinate system to define 
location ( = 0, 1 or 2), coordinates X-Y-Z (or R-0-2 or R-0-4>)/ and 
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Figure 4-1. INPUTS Input Deck Setup 



a) BASIC CARTESIAN X-Y-2 


Z 




c) SPHERICAL R-0-0 


Figure 4-2. Coordinate Systems 
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tdenHfication number of coordinafe sysfem to define displacemenfs 
( ~ 1, 2 or other). (Coordinate' Identification number 0 specifies the basic 

Cartesian system, 1 specifies the basic spherical system, and 2 specifies a 
special Cartesian system. Note: The displacement coordinate system I.D. is 
part of the BOPACE format, but is simply read and printed by INPUTB; also, 
for a 2-dimensional problem, the third coordinate is input as blank or zero. 

The basic coordinate systems are illustrated in Figure 4-2.) 

Format (215, 3F10. 0, 15) 

Blank card after last structural node. 

Thermal node location data. For each thermal node give node I.D. number, 
identification number of coordinate system to define location ( =-0, 1 or 2), 
coordinates X-Y-Z (or R-0-Z or R-0-i>). (For a 2-dimensional problem, the 
third coordinate is input as blank or zero.) User has .option of one or two 
thermal nodes per card. 

Format (2 (215, 3F10.0)) 

Blank card after last thermal node. 

Thermal node time-temperature data. 

a. Number of thermal times (’>2), initial default temperature. 

Format (MO, FI 0.0) 

b. Time-temperature vectors. For each vector; 

Time 

Format (FiO.O) 
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Thermal node temperatures. For each node give node l.D, and temperature. 
User has option of from one to four thermal nodes per card. 

Format (4(110, FI 0.0)) 

Blank card after last value of each vector. 

Note: If temperature is not given for any node, it is assigned a value by the 

program. For the first thermal time, this value is the default temperature 

in input Item 1-4. For later times, it is the value from the preceeding. 
time. 

I“5« Structural times. Give time value for each time at which interpolated structural 
temperature output is desired. 

Format (FlO.O) 

After last structural time, insert a card with an integer 9 in column 20. 

C-2 Blank card after last problem. (Multiple problem data sets are run consecutively, 
by stacking each data behind the pre'vious one). 
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5.0 SUMMARY. OF OUTPUT 


The description of [NPUTB output is here divided into two parts. The first covers output 
which is primarily an echo check of the input data and also includes some intermediate 
calculated results. This output is contained only on the first output file (UOUTl). 

The second part of the output consists of the final interpolated structural results to be 
used with BOPACE. This output is contained on both the first and second output files 
(UOUTl and UOUT2). 

5..T ECHO CHECK OF INPUT DATA, AND 'INTERMEDIATE RESULTS 

“ The first page of INPUTB output for a problem contains the 80-character 
problem I.D. title input as item 1-1, 

Structural Node Definitions - The info'mation given in input item 1-2 is printed. Values 
are the structural node number, node I.D., location coordinate system number 
(0 - basic Cartesian, 1 = basic cylindrical, 2 = basic spherical), the location coordinates 
(X-Y-Z, R-0-Z, or R-0-$), and the direction coordinate system number (0 = basic 
Cartesian, 1 = basic cylindrical, 2 =. basic spherical, >2 = I.D. of special system). 

Thermal Node Definitions - The information given in input item 1-3 is printed. Values 
are the thermal node number, node i.D., location coordi note system number (0 = basic 
Cartesian, 1 = basic cylindrical, 2 = basic spherical) and the location coordinates 
(X-Y-Z, R-0-Z, or R-0-«D). 
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Thermal Node Time ~ Tsmperat'ure Veci'ors ~ ~ The information given in input item 1-4 is 
printed. First the number of thermal times (temperature vectors) is printed, along with the 
initial default temperature. Then for each given thermal time is printed the time value, 

. followed by the list of node I.D.'s and associated temperature values for each specified 
node . 

Intermediate Interpolation Results - This output is provided for the user so that he may review, 
if desired, the thermal nodes and corresponding interpolation weighting factors which 
are associated with each structural node. Values printed are the structural node I.D., 
and associated four thermal node point I.D.'s and corresponding weighting factors. If 
only three thermal nodes were used to interpolate temperatures for a particular structural 
node, then the fourth weight is printed as zero and the fourth thermal node point is printed, 
equal to the first. Similarly if only two thermal nodes were used, the third weight is 
also printed as zero and the third thermal point is also printed equal to the first. If only 
one thermal point could be used, all weights but the first are printed as zero and all thermal 
points are printed equal to the first. 

5.2 FINAL RESULTS 

These are the final interpolated structural temperatures at each specified structural 
time, and are written onto both the first and second output files, in BOPACE format. 

For each structural time the time value Is printed on the first line or card (with the 
remainder of the card filled with asterisks to help the user identify the first card). 

Following this, is the vector of structural node I.D.'s and corresponding temperature 
values, for all structural nodes. At the end of each vector is printed a blank card as 
required for direct input to BOPACE. 
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6.0 SIZE LIMITATIONS 


ihe following variables are used to specify maximum size limitations in INPUTS. The 
values set for these variables are given In Table 6-1. 


NMAX2 
NMAX4 
NMAX52 = 
NAAAX54 = 
NMAX81 = 


maximum number of structural nodes, 
maximum structural node I.D. number, 
maximum number of thermal nodes, 
maximum thermal node I.D. number, 
maximum number of thermal node times. 


TABLE 6-1: INPUTS MAXIMUM SIZE LIMITATIONS 

VARIABLE VALUE 


NMAX2 

1,000 

NMAX4 

5,000 

NMAX52 

500 

NAMX54 

2,500 

NMAX81 

100 
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7.0 INPUTB ERROR MESSAGES 


INPUTB uses the FORTRAN STOP codes described In this section to Indicate error conditions 
which may occur during execution of the program. The following are explanations of the 
error STOP codes, listed by subroutine in which they occur. 

READO 

9999 Normal program exit (not an error) caused by reading final blank card after 

all problems are run. 


READM 

•701 

702 

704 

705 

READMT 

801 

802 

804 

805 

READT 

901 

902 

TCOMP 

1001 


Structural node I.D. exceeds maximum. 

I.D. of a structural node location coordinate system not equal to 0, 1 or 2, 
Number of structural nodes exceeds maximum. 

No structural nodes input. 

Thermal node I.D. exceeds maximum. 

I.D, of a thernal node location coordinate system not equal to 0, 1 or 2. 
Number of thermal nodes exceeds maximum. 

No thermal nodes input. 

Undefined thermal node I.D. used to specify temperature at a thermal time. 
Number of thermal times is less than 2 or exceeds maximum. 

Structural time is outside range of thermal times, or times are not In 
increasing order. 
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8.0 EXAMPLE PROBLEMS 


The example problems provided are infended I’D familiarize the user with the INPUTS 
program. Two problems are provided, a 2-dImefisIonal problem which introduces INPUTS 
in a more easily visualized 2-dimensional format, and a 3-dimensional problem which 
performs essentially the same functions in three dimensions. 
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8.1 


2-DIMENSlONAL PROBLEM 


The mesh used in this problem is shown in Figure 8.1-1 . The closed circles correspond 

to the structural nodes and^ the open circles to the thermal nodes. Temperature values, 

at TIME = 1.0, were assigned to the thermal nodes using the equation: 

T = 50x - 100/ 

A listing of temperature values is shown in Table 8.1-1 . 

Certain precautions were taken in arriving at input data for this example problem, in 
order to most effectively illustrate the interpolation procedure, and to guarantee that 
identical results would be obtained on different computer systems in the presence of 
small roundoff errors. For these reasons it was important that certain thermal points 
should not coincide exactly with any structural points, and that two thermal nodes and 
a structural node should not lie on the same line. It was also important that each 
structural node be enclosed within a triangle formed by three thermal nodes, which, 
although not always possible, increases the accuracy of the interpolation. Of course, 
the user running an actual problem needs to be less concerned with precautions such as 
those described here. 
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STRUCTURAL NODES 
THERMAL NODES 


Table 8.1-1. Thermal Node Temperatures for 2-D Problem at Time =1.0* 


Node 

Temperature 

Node 

Temperature 

1 

50 

20 

0 

2 

0 

9 

-250 

3 

125 

10 

-200 

4 

175 

8 

- 75 

5 

50 

11 

-150 

6 

-125 

40 

- 50 

7 

- 75 

12 

-325 

15 

-125 

13 

-300 

30 

75 

14 

-175 


*Note: Default Temperature at time 0.0 = 25,0. 
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START 

INPUTS 

5 

CHECKOUT 

6 7 

PROBLEM ( 2 -D 

MESH ) 


1 

0 0.000 

0.000 

0.000 


2 

0 1.000 

0.000 

0.000 


3 

C 2,000 

' 0.000 

0.000 


4 

0 3.000 

0.000 

0.000 


5 

0 d. 000 

l.OOO 

0.000 


6 

0 l.OOC 

1.000 

0.000 


7 

2 2.236 

26.525 

0.000 


8 

0 3.000 

1.000 

0.000 


9 

0 o.doc 

2,000 

0,000 


10 

0 1.000 

2.000 

0.000 


100 

0 2.000 

2.000 

0,000 


110 

0 3.000 

2 . 000 

0.000 


150 

0 0.000 

3,000 

0.000 


300 

0 1,000 

3.000 

0.000 


15 

0 2.000 

3 . 000 

0.000 


16 

• ■ 1 4.242 

45.000 

0.000 

00 

. 

1 

0 0.000 

- 0.500 

0.000 

1 

CO 

3 

0 1.500 

“ 0,500 

0.000 

5 

0 2.000 

0.500 

0.000 


7 

0 1.500 

1.500 

0.000 


15 

2 l .*581 

71.565 

0.000 


40 

0 4,000 

2.500 

0.000 


8 

0 2.500 

2 . 000 

0.000 


10 

0 1.000 

2.500 

0.000 


13 

0 2.000 

4.000 

0.000 


0.0 


2 0.0 


:.o 


.25 

.50 

1.0 


1 50.00 
6 ' - 125.00 
9 - 250.00 
40 - 50.00 

0 

0 

0 

9 


3 125.00 
7 - 75.00 
10 - 200.00 
12 - 325.00 


08 / 15/75 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


2 0 1.000 

4 0 3.500 

6 0 - 0;500 

30 0 3.500 

20 0 3.000 

12 0 -01 500 

9 0 0.000 

11 0 2.500 

14 1 4.950 


0.500 0.000 
0.000 0.000 
1.000 0.000 
1.000 0.000 

1.500 0.000 

3.500 0.000 

2.500 0.000 
2.750 0.000 
45.000 0.000 


4 175.00 
15 - 125.00 
8 - 75.00 
13 - 300.00 


5 50.00 
30 75.00 
11 - 150.00 
14 - 175.00 
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** NODE ** 

Nl7. —I. {).— LOCATE XJ -X?- X3 SPLAC 6 

1 1 0 0.0 0.0 0.0 0 

22 0 O.IOOOOE 01 0.0 0.0 0 

^ O 0.20000E- 01” 0,0 “ 0.0 0 

4 A 0 0.T00006 01 0.0 O.o 0 

55 0 0,0 O.lOOOOf 01 0.0 0 

— -0 #. 0 O.U'OOOt- 01”0.10000fc“01 -0.0 0 

7 7 ? 0.??3O0r 01 0,20'y^bL 02 0.0 0 

!) r. 0 o.iOOiiiJt 01 o.ioooot tn (uo , o 

— 4, — (, 0.0 - n.’ooooE <ii- 0.0 0 

10 10 0 0.100001 01 O.aOOOOE 01 0.0 0 

, 11 100 0 o.ioooor 01 o,;;(V4i((ul- oi o.o o 

, — , 12 . 110 -, 0 0 . 3 f '(.loot- 01 ~0.20(>{ifrt (‘1 -- 0.0 0 

I'i 1‘SO O l;.U 0.'}0(;ll0L 01 0.0 0 

lA ?,CiO (J G.lOoClit; 01 0.3OOUOF- 01 0.0 0 

in- . ir. 0 ().2 000C.!-“01 — 0.3()00()b- -01--0.0- 0 

10 10 1 0.m?a:-OF 01 O.Ai-OOoE 02 0.0 0 


I 

■-=N- 


■S' T NODE 

NO.- ■ I.D.- --LOCATE XI- — — ■ X2--- X3 

1 1 0' 0.0 -0.50000E 00 0.0 

2 2 0 n.l'onOOE 01' O.^iOOOOE 00 O/.O 

(f f). 1 501106 ■ 01-—-0.5OUOO& -yO - 0-0- 

i, /, 0 O.S.‘'OOOC ol 0.0 0.0 

'5 5 0 o.rooooE m o.^odooe no o.o 

^ ^ 0 O.EOOOoi 00 - 0. 1C0006- 01--0.0 

7 7 0 0.1''0ouE Ol 0.1500O6 01 0.0 

& 30' 0 O.ai'OOOE 01 O.lOOOOf (:1 0.0 

t, !<, 1 u-l*^fc;10E 01— Oi71‘6*>6 02 -0.0-™ 

10 20 0 O.jOuCjuE 01 O.lKfOoE ol 0.0 

11 40 0 0.400006 01 0.25000E 01 0.0 

12 12 0 0.500O0E- (HV--0.3500(»6 •ol — O.O- 

13 I' 0 o..?:.oont oi 0.20000F 01 0.0 

14 ^ 0 0.0 0.3‘-'00(,6 01 0.0 

15 10 fj O.lOOUOt-01 — 0.25('006 01 0.0 

Ih 11 0 0.2L0(u.>r Ol 0.17L00L Cl O.o 

17 IT 0 0.2i>iHniE 01 0. '•00006 01 O.U 

IK — 14 01 — ().'i5000h {>? 0,0 
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8.2 


3-DIMENSIONAL PROBLEM 


The strucf-ural network for this problem consists of three parallel planes, with the structural 
grids stacked one on top of the other, as-shown in Figure 8,2-1 . The arrangement of the 
nodes on each respective plane is shown in Figure 8,2-2. The closed circles correspond 
to the structural nodes and the open circles to the thermal nodes. The temperature values 
at TIME =1.0 were assigned to the thermal nodes using the equation: 

T = 50 X - lOOy + 50z 

A listing of the temperature values is shown in Table 8.2-1. 

The:same precautions were taken when formulating input data for this example problem as 
werei-aken for the 2-dimensional example, with additions being made fora 3-dimensional 
format . 

Thermal nodes were placed on levels .999, 1.999, and 3.001 so that thermal nodes and 
structural nodes would always be in different horizontal planes. These steps were utilized 
in order to avoid any inconsistencies in output data due to roundoff on different computer 
systems . 

It was also necessary to try to place the thermal nodes in such a way that each structural 
node v/as enclosed within a tetrahedron formed by four thermal nodes. When this condition 
could not be met the program would substitute the node found for point one in place of 
the missing points. This case is evident for nodes 12, 13, 30, etc., in the output data. 


8.2-1 







Figure 8.2-1. 3-D Problem Geomeuy Definition 


Table 8.2-1. Thermal Node Temperatures for 3-D Problem at Time =1.0* 


Node 

Temperature 

Node 

Temperature 

1 

24.95 

15 

124.95 

2 

25.05 

16 

-100.10 

3 

- 75.05 

17 

75.00 

4 

0.00 

18 

- 75.10 

5 

62.45 

19 

- 25.15 

6 

-124.95 

20 

-150.00 

7 

- 75.05 

21 

175.15 

8 

-175.15 

22 

75.05 

9 

- 75.15 

23 

25.05 

10 

25.50 

24 

- 50.00 

11 

174.95 

25 

50.10 

12 

175.10 

26 

-174.95 

13 

- 0.05 

27 

-124.95 

14 

25.0 

1 

• 


*Note: Default temperature at time 0.0 = 0.0. 
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8.2-4 


START 5" 6 7 


INPUTS 

10 
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0 
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1.000 
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0 
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1.000 
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20 

0 

l.OOO 

1 .000 
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1.000 
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s 
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9 
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PROGRAMMER MANUAL 



9.0 SUBROUTINES 


All 1-he subroutines of the INPUTB program are described in this section, with their 

calling sequence and argument variable definitions. The main program is described 
in Section 3. 

9.1 READO 

Subroutine READO reads the problem identification card, and the input/output file 
unit numbers. 


Call READO (UIN, UOUTl, UOUT2) 


UIN 

- 

input file unit number (e.g., cards). 

UOUTl 

- 

first output file unit number (e.g., printer). 

UOU.T2 

- 

second output file unit number (e.g., punch) 

9.2 

READM 


Subroutine READM reads the structural node locations, using BOPACE 3-D data format. 
Call READM (UIN, UOUTl, NOD, COORD, NODE, NMAX2, NMAX4) 


UIN 

input file unit number. 

UOUTl 

first output file unit number 

NOD 

number of structural nodes. 

COORD 

COORD (J, 1) = Jth coordinate of node 1. . 

NODE 

NODE (1) = external node l.D. for internal node I. 

NMAX2 

maximum number of nodes. 

NMAX4 - 

V 

maximum node l.D. number. 
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9.3 


READMT 


Subroutine READMT reads the thermal node locations. 

Call READMT (UIN, UOUTl, NODT, COORDT, NODET, NODIT, NMAX52, NMAX54). 
UIN - input file unit number. 

UOUTl - first output file unit number. 

NODT - number of thermal nodes. 

COORDT ” COORDT (J, I) =J th coordinate of node I . 

NODET - NODET (I) = external node I.D. for internal node I. 

NOD IT “ NODIT (I) = node internal nurrtier for external I.D. I. 

NMAX52 - maximum number of thermal nodes. 

NMAX54 - maximum thermal node I.D. number. 


9.4 READT 

Subroutine READT reads the vector of thermal node temperatures, at each given thermal 
time. 


Cal! READT (UIN, UOUTl, lOT, NTIME, TIME, NODT, NODIT, NMAX54, NMAX81) 


UIN 

UOUTl 

lOT 

NTIME 

TIME 

NODT 

NODIT 

NMAX54 - 

NMAXOl '■'i.'T 


input file unit number 
first output file unit number. 

file unit number for storing thermal node temperature vectors, 
number of temperature vectors (time values) given for thermal nodes. 
TIME (I) = I tht ime value for thermal nodes, 
number of thermal nodes. 

NODIT (I) - node internal number for external I.D. I. 
maximum thermal node I.D. number. 

maximum. number bf thernial node temperature vectors,(tjme Yffliyes),,,, 
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9.5 


COMPW 


Subroutine COMPW selects four thermal. nodes and associoted weighting factors for each 
structural node, and stores them for later interpolation. ’ 

Call COMPW (UOUTl, NOD, NODT, COORD, COORDT, NON, NOW, NODE, NODE!) 
UOUTl - first output file unit number 

nod - number of structural nodes, 

NODT - number of thermal nodes, 

COORD “ COORD (J, I) - Jth coordinate of structural node I. 

COORDT - COORDT (J, 1) = Jth coordinate of thermal node 1. 

non - NON (J, I) = selected jth thermal node (J= 1-4) for weighting 

temperature of structural node 1. 

NOW - NOW (J, I) = computed J the thermal node weight (J = 1-4) for 

structural node I . 

NODE - NODE (I) = external l.D. for internal structural node }. 

NODET - NODET (I) — externa! l.D. for internal thermal node I. 

9.6 TCOMP 

Subroutine TCOMP perfcrms the actual interpolation in space and time, to compute 
structural node temperatures, using the thermal node temperature vectors and the 
weighting factors. 

Call TCOMP (UIN, UOUTl, UOUT2, lOT, NOD, NODT, NON, NOW, NTIME, 

TIME, NODE) 

UIN - input file unit number 

UOUTl - first output file unit number. 

UOUT2 - second output file unit number. 
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[OT 

NOD 

NODT 

NON 

NOW 
NT I ME 

TIME 

•NODE 


file unit number for storing thermal node temperature vectors, 
number of structural nodes, 
number of thermal nodes. 

NON (J, I) - Jth thermal node (J = 1-4) for .weighting temperature 
of structural node I. 

NOW (J, I) = J the thermal node weight (J = 1-4) for structural node 1. 
number of time values (temperature vectors)- for thermal nodes. 

TIME (I) = I th time value for thermal nodes* 

NODE (I) = external I.D. for internal structural node I. 
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10.0 COMMON BLOCKS 


The only common used in the INPUTB program are the two labeled common blocks 
COMTO and COMTl, which provide storage space for thermal node temperature vectors. 

Common/COMTO/TO (used in MAIN, TCOMP) 

Common/COMTl/Tl (used in MAIN, READT, TCOMP) 

The MAIN program contains both TO and Tl, and is used to dimension these vectors. 
Subroutine READT uses the Tl vector as temporary storage to read in each thermal node 
temperature vector before writing it onto the file lOT. Subroutine TCOMP uses both 

the TO and Tl vectors, to read and store thermal node temperature vectors from file 

lOT, at successive thermal time values. These vectors are then used to interpolate all 
structural node temperatures within the time interval from TO to Tl. 
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11.0 FILE USAGE 


INPUTB uses FORTRAN I/O fo access several files. Some of the files are defined by the 
user and the others are defined In the program. A list of files by logical unit number 


fol 1 ows: 



UNIT NUMBER 

DESCRIPTION 

DEFINED BY 

5 

Input card file 

READO 

UIN 

Input data file 

user 

UOUTl 

First output file. Contains all output 
results, including echo check of Input 
data, and intermediate results. 

user 

UOUT2 

Second output file. Contains only 
final output results, i.e., interpolated 
structural node temperature data. 

user 

lOT (=11) 

Storage file for thermal node temperature 
vectors. 

MAIN 
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APPENDIX - INPUT B PROGRAM LISTING 





I N P U r B (THERMAL/STRUCTURAL DATA INTERFACE PROGRAM) 00000020 
INPUTS IBM 360/370 VERSION (1000 NODES) DATED 09/15/75 00000030 
BOEING AEROSPACE COMPANY, P.O. BOX 3999, SEATTLE, WASHi 98124 OOOCOOAO 
ENGINEER/PROGRAMMER: R*G. VOS, PHONE 773-2946, BLOG. 18-05 00000050 

♦ :<c » 5j< « 4; <5 :(J :{, >}c ;^c T{c ««>)(** 0 n 0 0 0 0 6 0 


SPACE-TIME INTERPOLATION FOR 1-, 2- OR 3-DIMENSIONAL SPACES* 
PROBLEM SIZE CAPABILITY INCLUDES 1000 STRUCTURAL NODES, 

500 THERMAL NODES, 100 THERMAL NODE TIMES, ARBITRARY NUMBER OF 
STRUCTURAL NODE TIMES. 

MAXIMUM STRUC rURAL, THERMAL NODE I.D. NUMBERS ARE 5000,2500* 
NOCES ARE LOCATED IN BASIC CARTE S I AN , C YL I NORI CAL OR SPHERICAL 
COORDINATES. 

integer UIN,UDUT1,U0UT2, 10T,NOD,NODT,NTIME,NMAX2,NMAX4,NMAX52, 
INMAX5A, NMAX8I 

integer NCDEUOOO),NOOET(500),NODIT{250D>,NON(4, 1000) 

REAL T0( 500) ,T1 ( 500) , T IME ( 100 ) , C00RDI3 , 1000 ) ,COORDTi 3,500 ) , 
1NOW(4,1000) 

COMMON/COf'TO/TO/COMTl/Tl 
DATA lOT/ll/ 

DATA NMAX2,NMAX4,NMAX52,NMAX54, NMAX81 / 1 000 , 5000, 500 ,2500 , 100/ 
READ PROBLEM I.D. AND FILE UNIT NUMBERSi 
. CALL REAOO(UIN,UOUT1,UOUT2) 

READ STRUCTURAL NODE DEFINITIONS. 

CALL READM(UIN,UOUT 1,N0D, COORD ,NODE , NMA X2 , NMA X4) 

READ THERMAL NODE DE F I N I T TONS . 

CALL READMT(UIN,U0UT1 ,NDDT,C00RDT,N0DET,N0DIJ,NMAX52,NMAX543 
READ THERMAL NODE TIME-TEMPERATURE VECTORS. 

CALLREADT( UIN,UQUT1 , 10 T, NT IME, T I ME , NODT , NOD I T, NM AX54 , NMAX8 1 ) 


C0OCOO7O 
00000080 
00000090 
COOCOlOO 
00000110 
00000120 
C00C0130 
OOOOOIAO 
00000150 
00000160 
00000170 
00000180 
00000190 
00000200 . 
00000210 
00000220 
00000230 
00000240 
00000250 
00000260 
00000270 
00000280 
00000290 
00000300 


SELECT THERMAL POINTS AND COMPUTE CORRESPONDING WEIGHTS TO GIVE 00000300 
STRUCTURAL TEMPERATURES AS FUNCTIONS OF THERMAL NODE TEMPERATURES. 00000310 
CALL COMPWtUOUTl , NOD, NODT, COORD, COOROT, NON , NOW, NODE, NODE T) 00000320 

COMPUTE AND OUTPUT STRUCTURAL NODE TEMPERATURES AT GIVEN TIMES. 00000330 
CALL TC0MP(UIN,U0UT1,UDUT2, IOT,NOD,NODT,NON,NOW,NTIME,TIME,NOOE> 00000340 
^ 00000350 

END 


c 


SUBROUTINE RE ADO ( U I N ,U0UT1 ,U0UT2 ) 

READ FILE UNIT NUMBERS AND PROBLEM IDENTIFICATION 


00000370 

00000380 



3mGINAL PAGS'H 
OF POOR QUAUT3 


C 

C 

C 


UIN = INPUT FILE UNIT NUMBER. 

UOUTltUOUT2 = OUTPUT FILE UNIT NUMBERS. 

STOP 9999 = NORMAL PROGRAM EXIT AFTER LAST PROBLEM IS INPUT. 
integer UIN,U0UT1,U0UT2 
INTEGER I ,START,STAR, IDENT (20) 

101 rORMAT(AA,6X,3I5) 

102 FGRMAl (20AA ) 

2Cl FtlRMATI 1H1,20A4) 

DATA STAR/4HSTAR / 

REACISt 10USTART,UIN,U0UT1,U0UT2 
IF(STARr.NE.STAR)STOP 9999 
REAOIUIN, 102) (IDENTII ) tl=lt20) 

IFIUOUTl .GT.0)WRITE(U0UT1t 201 ) { IDENTI 13,1 = 1,2 0) 
IF(UCUT2.GT.0) WRITE(U0UT2r 102) I I DENT I I ),.I=1,20) 

RETURN 

END 


00000390 

00000400 

00000410 

00000420 

00000430 

00000440 

00000450 

000.00460 

00000470 

00000480 

00000490 

00000500 

00000510 

C0000520 

00000530 

00000540 




SUBROUTINE RE ADM I UI , UO , NOD , COORD ,N00£ ,NMAX2 , NMAX4 ) 

C READ STRUCTURAL NODE DATA.' 

C UI,UO = INPUT, OUTPUT FILE UNIT NUMBERS. 

C NOD = NUMBER OF STRUCTURAL NODES, 

C CGORDIJ,!) = COORDINATES OF NODE I. 

C NODE! I) = NODE EXTERNAL I.D. FOR INTERNAL NUMBER I; 

C NMAX2 = MAXIMUM NUMBER OF NODES. 

C NMAX4 = MAXIMUM NODE I.O. NUMBER. 

C STOP 701 = STRUCTURAL NODE I.D. EXCEEDS MAXIMUM. 

C STOP 702 = I.D. OF A STRUCTURAL NODE LOCATION COORDINATE SYSTEM 

C NOT EQUAL TO 0, 1 OR 2. 

C STOP 704 = NUMBER OF STRUCTURAL NODES EXCEEDS MAXIMUM. 

C STOP 705 = NO STRUCTURAL NODES INPUT, 

INTEGER UI, UO, MOD, NODE! 1) ,NMAX2,NMAX4 
REAL C00RC(3,1) 

INTEGER I ,LC00RD,DC00RD 
REAL F, ANGLE, X,Y, Z 
101 FORMAT (2I5,3F10.0, 15 ) 

201 FORMAT! 1H1,10H«# NUDE *=<'/lH ,18H NO. liD. LOCATE, 6X, 

12HX I, IIX, 2HX2, IIX, 2HX3, 8X, 8HDISPLACE) 

202 FORMAT (2 15, 2X, 15, 3X , 3 ( E12 .5, IX ) , IX, L5) 

F = 3.1415927/180. 


00000550 

0C000560 

00000570 

00000580 

00000590 

00000600 

00000610 

00000620 

00000630 

00000640 

00000650 

00000660 

00000670 

00000680 

00000690 

00000700 

00000710 

COOC072O 

00000730 

00000740 

00C00750 

00000760 



ooooooooo 


WRITE(UO, 201 ) 

NCD = 0 

6 READJUI, 1‘Cl) I »LCOORD,X,YfZ,DCOORO 
IFd.LEiOJGO TO 150 
IF( I.GT.NYAX4)ST0P 701 

IF ( LCOORO.NE. 0. AND. LCOORD.NE. 1 . AN0.LC00RD.NE.2 )ST0P 702 
NOD = NOD+1 

IF(NGD.GT.NMAX2)ST0P 704 

IF{ UO.GT.O} WRITE(U0,202 )N0D, I,LCOORDtX,Y,Z,DCOORD 
NODE (NOD) = I 
IFaCCORD.EQ.OGO TO 12 
IF( LCOORD.EQ. 1 )GO TO 11 
C SPHERICAL COORDINATES. 

ANGLE = Z*F 
Z = X*SINIANGLE) 

X = XX^CUS(ANGLE) 

C CYLINDRICAL COORDINATES. 

II ANGLE = Y«F 

Y = X^C'SINIANGLE) 

X = X=4'C0S I ANGLE) 

C BASIC CARTESIAN COORDINATES. 

12 CQORDt I, NOD) = X 
C00RD(2,N0D) = Y 
CCORDO.NCD) = Z 
GO TO 6 

150 IF(NOD.EQ.O)STOP 705 
RETURN 
END 


00000770 

CCC00780 

00000790 

00000800 

00000810 

00000820 

00000830 

CC000840 

00000850 

00000860 

00000870 

00000680 

00000890 

00000900 

00000910 

00000920 

00000930 

000C0940 

00000950 

00000960 

00000970 

00000980 

00000990 

00001000 

00001010 

00001020 

00001030 

00001040 


SUBROUTINE READMTIUI » UO, NODT , COORDT , NODET » NOD IT, NMAX52, NMAX54 1 
READ THERMAL NODE DATA. 

UI,UO = INPUT, OUTPUT FILE UNIT NUMBERS. 

NCDT = NUMBER OF THERMAL NODES. 

COORDTIJ,!) = COORDINATES OF NODE I, 

NOOET(I) = NODE EXTERNAL I.D. FOR INTERNAL NUMBER I. 

NCDIT(I) = NODE INTERNAL NUMBER FOR EXTERNAL L.D. I. 

NMAX52 = MAXIMUM NUMBER OF THERMAL NODES. 

NMAX54 = MAXIMUM THERMAL NODE I.D. NUMBER. 

STOP 801 = THERMAL NODE I.D. EXCEEDS MAXIMUM. 


00001050 
00001060 
00001070 
C0001080 
00001090 
OOOOllOO 
OCOOlllO 
00001120 
00001130 
00001 140 
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STOP 802 = I.D. of' a THERMAL MODE LOCATION COORDINATE SYSTEM 
NCT EQUAL TO 0,1 DR 2* 

STOP 804 = NUMBER OF THERMAL NODES EXCEEDS MAXIMUM. 

STOP 805 = NO THERMAL NODES INPUT. 

INTEGER UI,UO,NODT,NODET( I ) ,NODITI 1 )',NMA X52 , NMA X54 
REAL C00RDT(3,1) 

INTEGER I ,LCOCRO,K, IST0R(2 ) ,LST0R(2) 

REAL Ft ANGLEfX, Y,Z,XSTOR< 2) ,YST0R( 2) ,ZST0R( 2) 

101 FORMAT (2 12I5,3F10.0) ) 

2Gl FORMAT! IHl, 10H<‘ T NODE */lH ,18H NO. UD. LOCATE, 6X, 

12HX1, 11X,2HX2, 1IX,.2HX3) 

202 F0RMAT(2I5,2Xt I5,3X , 3 ( E12.'5, IX) ) 

F = 3.1415927/180. 

W.RITE(U0,201) 

DO 2 I=l,NMAX54 
2 NOD IT ( U = 0 
NODT = 0 

6 READ (Ul, 101 ) { ISTORIK),LSTOR(K ) ,XSTOR(K) , YSTORIK) ,ZSTOR(K) ,K=l,2) 
DC 8 K-=l,2 

IFU STOR(K) .NE.OGO TO 9 

8 CONTINUE 

IF{NCDT;EG.0)ST0P 805 
RETURN 

9 DO 20 K=l,2 

• I = ISTOR(K) 

IF{ L.LEi.0)G0 TO 20 
IF( I .GT.NMAX54)ST0P 801 
LCCCRD = LSTQR(K) 

IF ( LCOORD.NE.O. AND. LCQORD . NE . 1 . AND. LCOORD .NE.2 ) STOP 802 
NODT = NODT+1 
IF(N0DT.6T.NMAX52 )STOP 804 
X = XSTOR(K) 

Y = YSTORIK) 

Z = ZSTOR(K) 

IFIUO.GT.O) WRITE (UO, 202 ) NODT, L, LCOORD, X„Y,Z 

NOOET(NOOT) = I 

NODIT ( I ) = NDDT 

IF{LCOORD.EO.O)GO TO 12 

IF (LCOORD. EQ. UGO TO 11 

SPHERICAL COORDINATES. 


00001150 
00001160 
C0001170 
00001180 
00001190 
00001200 
00001210 
00001220 . 
C0001230 
00001240 
00001250 
CC001260 
00001270 
00001280 
00001290 
00001300 
00001310 
00001320 
00001330 
00001340. 
00001350 
00001360 
00001370 
00001380 
00001390 
00001400 
00001410 
00001420 
0000] 430 
00001440 
CC001450 
00001460 
00001470 
00001480 
00001490 
00001500 
00001510 
CC001520 
00001530 
00001540 



nnoooooooono 


ANGLE = Z«F 
I = X’f'S IN (ANGLE) 

X = X*COS (ANGLE) 

C CYLINDRICAL COORDINATES- 

11 ANGLE = Y*F 

Y = X*SIN(ANGLE) 

X = X*COS(ANGLE) 

C BASIC CARTESIAN COORDINATES. 

12 COORDT(1,NODT) = X 
CCGRDT(2tNOOT) = Y 
COORDTI 3,NODT) = Z 
20 CONTINUE 
GO TO 6 
END 


C0001550 

00001560 

00001570 

00001580 

00001590 

00001600 

00001610 

00001620 

00001630 

OC001640 

00001650 

00001660 

00001670 

OOOOL680 


101 

102 

103 

201 


SUBROUTINE READTIUI »U0» lOTtNT IME,T IME,NODT»NODITt NMAX54,NMAX81) 
READ THERMAL NODE TIME-TEMPERATURE VECTORS. 

UltUO = INPUT, OUTPUT FILE UNIT NUMBERS. 

lOT = FILE UNIT NUMBER FDR THERMAL NODE TIME-TEMPERATURE VECTORS. 
NTIME = NUMBER OF TIME VALUES FOR THERMAL NODES; 

TlMEd) = ITH TIME VALUE FOR THERMAL NODES. 

NCDT = NUMBER OF THERMAL NODES. 

NODIT(I). = THERMAL NODE INTERNAL NUMBER FOR EXTERNAL UD. I.' 
NMAX54 = MAXIMUM THERMAL NODE L.D. NUMBER. 

NMAX81 = MAXIMUM NUMBER OF THERMAL TIMES. 


STOP 901 = UNDEFINED THERMAL 
STOP 902 - NUMBER UF THERMAL 
MAXIMUM. 

COMMON/ CQMTl/Tl 

TKl) = TEMPERATURE OF ITH THERMAL NODE AT GIVEN TIME. 
integer UI,U0,U0T, NTIME, N0DT,N0DIT( 1 ) , NM AX54, NMA X 8 1 
REAL TIME(l) ,T1( 1) 

INTEGER I,K,ITIME,ISTOR(A) 

REAL TEMPO, ST0R(4) 

FORMATI IlOtFlO.O) 

FORMAT(FIC.O) 

F0RMAT(4{ 110, FIO.O) ) 

FCRMAT(1H1,40HNUMBER OF TIMES (TEMPERATURE VECTORS) = ■ 
llCX,22HnEFAULT TEMPERATURE = ,F10.4) 


00001690 

00001700 

00001710 

00001720 

00001730 

00001740 

00001750 

00001760 

00001770 

C0001780 


NODE I.D. USED TO SPECIFY TEMPERATUREOOOOl 790 
TIMES IS LESS THAN 2 OR EXCEEDS 00001800 

00001810 
00001820 
00001830 
CC001840 
00001850 
00001860 
00001870 
00001880 
00001890 
00001900 
13, 00001910 

00001920 






202 F0RI^ArUHl»30HTEMPERArURE VECTOR FOR TIME = ,E12.A/ 
11H0,25HTHERHAL PT. TEMPERATURE) 

203 FORMATdH 1 1 5 , 7X , F 12 . 2 ) 

REAC(UI, 101 )iMTIME,TEMPO 
IF(UO.GT.O) WRITE (UO, 201 )NTIME,TEMPO 

IF ( NT IME.LT.2.0R.NTIME .GT.'NMAXSl) STOP 902 
REWIND lOT 

DC 100 ITIME=l,NTIME 
READ(UI,.1C2)TIME( ITIME) ' 

IF (UO.GT.O)WRITE(UO, 202 )TIME( I.TIME) 

IF (IT I ME. GT .1 ) GO TO 13 
DO 5 I=1,N0CT 
5 Tim = TEMPO 

13 READ (UI ,103) ( I STOR (K ) , STOR ( K ) , K= 1 ,.A ) 

DO 15 K=1,A 

IF( ISTOR(K) .NE.O)GO TO 16 

15 CONTINUE 
GO TO 99 

16 DO 20 K=l,4 

I = ISTOR(K) 

IF! L.LE.OGO TO 20 
IF(UO.GT.O)WRITE(UO,203) I ,STORIK) 

IF( I.GT.NMAX5A)ST0P 901 
I = NODIT(I) 

IF! I-LE.OSTOP 901 
Tim = STORIK ) 

20 CONTINUE 
GO rO 13 

99 WRITEUOn (Tim, I=l,NOOT) 

100 CONTINUE 
RETURN 
END 


C0001930 

000019AO 

000,01950 

00001960 

00001970 

00001980 

00001990 

CC002000 

00002010 

00002020 

00002030 

000020AO 

00002050 

C00C2060 

00002070 

00002080 

00002090 

00002100 

00002110 

00002120 

00002130 

00O021A0 

00002150 

00002160 

00002170 

oooo 2 iao 

00002190 

00002200 

00002210 

00002220 

00002230 

000022A0 


C 

C 

C 

C 

C 


SUBROUTINE COMPW ( UO , NOD, NODT, COORD tCOORDT, NON , NOW , NODE ,NODET) 
SELECT THERMAL POINTS AND SET WEIGHTS FOR EACH STRUCTURAL NODE. 
UO = OUTPUT FILE UNIT NUMBER. 

NOD = NUMBER OF STRUCTURAL NODES. 

NCDT = NUMBER OF THERMAL NODES. 

CGORD(J,I) = JTH COORDINATE OF STRUCTURAL NODE I. 


00002250 

00002260 

00002270 

00002280 

00002290 

00002300 



o r» o o o 


STRUCTURAL 

STRUCTURAL 


NODE I 
NODE I 




COOROT(J,.I) = JTH COORDINATE OF THERMAL NODE I.. 

NCN(J,I) = JTH THERMAL NODE NUMBER (J=l-4) FOR 
NOW(J,I) = JTH THERMAL NODE WEIGHT (J=l-4) FOR 
NODE! I) = EXTERNAL I.D. FOR STRUCTURAL NODE I. 

NCDET(I) = EXTERNAL I.D, FOR THERMAL NODE I. 

INTEGER UO»NOO,NODT »N0N(4» 1 ) ♦NODE! 1 ) ♦NDDETI I ) 

REAL COORD! 3, 1) ,COOROT( 3, 1) tNOWI^,!) 

INTEGER INOD, I,Jt I 1, I2f I3» 14 

REAL C,R,W,W1 ,W2,W3,W4»P(3),P0(3I,PU3)»P2( 3),P313),P4(3), 
lS(3)tS12(3),S23(3) 

201 F0RMATI1H1,5H N00E»1X,4(4H PT , 1 2 ) » 4( 6X, 2HWT , 1 21 ) 

2 02 FORMAT! 15, 2X, 4 ( 1 5 , 1 X ) , 2X , 4F1 0 .4 ) 

IF!UO.GT.OJWRITE!UO,201) I J,J=1,4) ,! J,J=1,4) 

00 1000 IN0D=1,N0D 
DC 5 J=l,3 

5 PO!J) = COORD! J , INOD) 

LOCATE 1ST THERMAL PT AS NEAREST PT 
R -= 10.E30 
DO 500 I=l,NODT 
DO 456 J=l,3 

456 P!J) = COCRDTIJ,!) - POIJ) 

IF! P{ i)>!'*2 + P!2)«=«'2 + P!3)**2.GE.R)G0 TO 500 

11 = I 

R = P ! 1 )**2+P ! 2 )*«2-*;P! 3 ) 4*2 
500 CONTINUE 

IFIR.EQ.O. )G0 TO 951 
DO 512 J=l,3 

512 P1!J) = CC0RDT!J,I1) - POIJ) 

LOCATE 2ND THERMAL PT AS NEAREST PT AT LEAST 90 DEGREES FROM 1ST 
R = 10.E30 

12 = 0 

DO 600 I=1,N00T 
DO 556 J=l,3 

556 P!J) = COCRDTIJ,!) - PO!J) 

IF! PI 1)4*2+P! 2)**2 + P!3)**2.GE.'R.0R.I.EQi IlOGO TO 600 
IF(P!1)4P1(1)+P(2)*PI!2)+P(3)*P1!3).GT.O,)GO TO 600 
12 = I 

R = PI 1 )**2+P(2)**2 + P!3)**2 
600 CONTINUE 

IF! I2.EQ.0)G0 TO 951 


00002310 
0C002320 
00002330 
00002340 
OC002350 
00002360 
00002370 
00002380 
00002390 
00002400 
00002410 
00002420 
00002430 
00002440 
00002450 
00002460 
00002470 
00002480 
00002490 
00002500 • 
C0002510 
00002520 
00002530 
00002540 
00002550 
00002560 
00002570 
00002580 
00002590 
00002600 
C0002610 
00002620 
00002630 
00002640 
00002650 
00002660 
00002670 
00002680 
00002690 
CC002700 



n o 


CO 


DO 612 J=l,3 

P2(J1 = C00R0T(J,I2) - PO(J) 

612 Si2 ( J) = P2 ( J ) - PI ( J ) 

LOCATE 3RD THERMAL PT AS NEAREST PT FORMING ENCLOSING TRIANGLE 
FOR PERPENDICULAR TO PLANE 
R = 10.E30 
13 = 0 

DO 700 I=1,N00T 
DO 656 J=It3 

656 P{J) = COORDT(J,I) - PO(J) 

IF(P( 1 )«*2 + P{ 2)**2+P( 3)*=«'2.GE..R.0R.I.EQ. I l.OR.:l.EQi 12)G0 TO 700 
DO 662 J=l,3 

662 S23(J) = P(J) - P2(J) 

S(l) = S12(2)*S23(3)-SI2{3)*S23(2) 

S(2) = S12(3 )«S23( 1 )-S12t 1 )«S23( 3) 

S(3) = S12m«S23(2)-SL2{2)*S23{l) 

IF(S( I).EQ.0,.AND.S{2)-EQ.0..AND.S(3).EQ.0,) GO TO 700 
W1 = S<1)«(P2(2)*P( 3)-P2( 3)*P( 2) ) tS( 2)«IP2(3)«PU)-P2( l)«p(3) ) 
l + S{3)*(P2m*PI2)-P2(2)*P<l ) ) 

W2 = S{ 1)#{P(2)*P1(3)-P(3)*P1(2) )+S(2)*{P{3I«P1(1)-P(1)#P1{3)) 
l+S(3J=5‘{Pm*Pl(2)-P(2)*Pl( 1) ) ^ uuuu^viu 

W3 = S ( 1 tPl {2)«P2 { 3 )-Pl (3 )«P2'( 2) ) +S ( 2 ) *( P 1 { 3 ) *P 2( 1 ) -P 1< 1)>5‘P2{ 3) ) 000 02 920 
l + S{3)*(Pim*P2(2)-Pl (2 )>i=P2(l ) ) 00002930 

IF((Wl.'LE.O.’«AND.W2.LE.O».'AND.W3«LE»0«)iOR. 00002940 

1 { W1 « G E • 0 . • AND • W2 • G E «0 • « AN D . W 3 . GE • 0 • ) ) GO TO 690 00002950 

00002960 

^ 00002970 

R = P(l)**2+P{2)«*2+P<3)*^2 00002980 

700 CONTINUE 00002990 

IF( 13.EQ.C)G0 TO 961 00003000 

DC 712 J=l,3 00003010 

712 P3(J) = CC 0RDT(J»I3)— PO(J) 0C003020 

C LOCATE 4TH THERMAL PT AS NEAREST PT FORMING ENCLOSING TETRAHEDRON 00003030 

TA~ 00003040 

It , 00003050 

DO 800 I=1,N0DT 00003060 

DO 756 J=l»3 00003070 

756 P(J) = COCRDT(Jfl) — P0(J) 0C003080 

IF(P( l)+«2+P(2)«>i‘2+P(3)>»=«2.GE.R.0R.'I.EQiI1.0R.I.EQ.I2,0R.l.EQ^I3) 00003090 
IGO TO eoo 00003100 


00002710 

00002720 

C0002730 

00002740 

00002750 

00002760 

00002770 

00002780 

00002790 

00002800 

00002810 

00002820 

00002030 

00002840 

00002850 

C0002860 

00002870 

00002880 

00002890 

00002900 

00002910 
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C 


W1 = P ( I ) * (P2(2 ) *P3 ( 3 )-P2 (3 ) *P3 ( 2 ) ) 

UP( 2)#(P2(3)*P3( 1)-P2( 1 )*P3( 3 ) ) + P { 3 ) « < P2 ( 1 ) «P3 ( 2 )-P2 ( 2 ) *P3 ( 1 1 ) 

W2 =~P1 ( 1 P3{ 2) *P ( 3 )-P3( 3 )*P J 2) ) 

1-P1(2)*(P3(3)*P11)-P3(1 )*P{3 ) )-Pi(3 )=<'(P3( l)«P(2)-P3l 2)*P( 1) ) 

W3 «= P2{ 1 )*( P( 2) *P1 ( 3) -P ( 3 )*P1 (2) ) 

1+P2 (2)*(P(3)*PUl>-P(l)^Pl(3))+P2(3)»{P(l)«Pl(2)-P(2)’i'PUl}) 

W4 =-P3 ( 1 ) «( PI I 2 ) «P2 ( 3 )-Pl ( 3 ) *P2 ( 2 ) ) 

1-P3(2)«( P1(3)*P2( D-PH 1 )>:‘P2(3) ) -P3 { 3 ) *{ P 1 1 1) *P2 ( 2 )~P H 2 ) #P2 ( 1) ) 
IF( ( Wl.LE.O..AMO.W2.LE.O.'.AND.W3.LE.O-. .AND.W4.L6.0.3.0R. 

I ( WL . GE .0 « - AND .W2 .G£ .0 w . AND.W3.GE.0 . . AND.W4.GE . 0. ) ) GO TO 790 
GO TO 800 
790 lA = I 

R = P ( 1 ) «=<<2+P ( 2 }«#2+P ( 3 ) 

8C0 CONTINUE 

IF{ IA.EQ.03G0 TO 971 
00 812 J=l,3 

812 P4U) = CC0RDTU,I4) - PON) 

4 THERMAL PTS CAN BE USED 

W1 - P4(1)*(P2{2)=«‘P3( 3)-P2t 3)«P3(2) ) 

1 + P4 (2)*( P2 13 )*P3 1 D-P2 ( 1 )«P31 3 ) )+P4( 3)*1 P2( 1 ) >4P3{ 2)-P2< 2)*P3( 1 ) ) 
W2 =-PU 1)*(P3(2)*P4(3)-P3(3)*P4(2)) 
l-Pl ( 2 )«1P3 ( 3)4P4( 1)-P3( 1 )«P4( 3) )-Pll 3)*1 P31 1 ) =«'P4 { 2) -P3 { 2 ) *P4( 1 ) ) 
W3 = P2 ( 1 )* ( P4 (2 ) =«'P1 1 3 )-P4 I 3 ) *P1 ( 2 ) ) 

1 + P2(2)*1P413)«P1(1)-P4(1)«P1(3))+P2(3)*(P4(1)*PU2)~P4(2)’<'P1{1)) 
W4 =-P3(l)*{Pl(2)*P2(3)-PH3)*P2(2)) 

l-P3{2)*(Pl(3)*P2 11}-PUl)»P2(3))-P3(3)#(Pltl)«P2(2)~PU2)*P2(l)) 

W = W1+.W2 + W3 + U4 

IFIW.EQ.O. )G0 TO 971 

W1 ='W1/W 

W2 “ W2/W 

W3 =: W3/W 

W4 = W4/W 

GO TO 991 

ONLY 1 THERMAL PT CAN BE USED 
951 12 = 11 
13 = II 
14-11 
W1 = 1.0 
W2 = 0. 

W3 = 0. 


00003110 

00003120 

00003130 

00003140 

00003150 

00003160 

C0003170 

00003180 

00003190 

0C003200 

00003210 

00003220 

00003230 

00003240 

00003250 

00003260 

C0003270 

00003280 

00003290 

C0003300 

00003310 

00003320 

CC003330 

00003340 

00003350 

000C3360 

00003370 

00003380 

00003390 

00003400 

00003410 

00003420 

00003430 

0'0003440 

00003450 

00003460 

00003470 

00003480 

00003490 

00003500 
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C 


W4 = 0. 

GC TO 991 

ONLY 2 THERMAL PTS CAN BE USED 
961 13 II 
14 = II 
W3 = 0, 

W4 = 0. 

W = 0. 

R = 0. 

DO 965 J=l,3 

W = W +' P2 (J }*(P2 (J )-Pl ( J ),> 

965 R = R -h IP2{ JJ-PK 
IFIR.EQiO. )G0 TO 951 
W1 = W/R 
W2 = 1.0-1>1 
GO TO 991 

ONLY 3 THERMAL PTS CAN BE USED 
971 14 = II 
W4 = 0. 

DO 975 J=l,3 

975 S23U) = P3U) - P2 (J ) 

S(l) = S12I2)*S23(3)-S12(3J*S23(2) 

S(2) = S12(3)<‘S23(1)-S12{ 1]«S23(3) 

S(3) = S12a)*S23(2)-S12(2)*S23(l) 

• IFISI n-EQ.0..AND.S(2).EQ.'0..AND,S(3).EQ.0.3 GO TO 961 


00003510 

CC003520 

00003530 

00003540 

G0003550 

00003560 

00003570 

C0003580 

00003590 

00003600 

0C003610 

00003620 

00003630 

00003640 

00003650 

00003660 

00003670 

00003680 

00003690 

00003700. 

00003710 

00003720 

00003730 

00003740 

00003750 


W1 = S( 1 )4( P2(2)*P3I 3)-P2(3)>;‘P3( 2) )+;SI 2)4(P2< 3)»«‘P3( 1 ) -P2 ( 1 ) *P3 1 3) ] 00003 760 
1+SI3)>!'(P2(1)>!‘P312)-P2(2)*P311 I ) 0C003770 

W2 = S( IH‘(P3(2)*P1(3)-P3(3)*P1{2) ) +S ( 2 ) * [ P3 ( 3 ) *P 1 1 1 J-P3 < 1) «P1 ( 3 )) 00003780 


l+S(3)=>^IP3(l)*Pl(2)-P3{2)«Pl{in 00003790 

W3 = S(1 )=frIPl (2 )«P2t 3 )-Pl (3 )#P2( 2) )+S{2)*(Pl( 3)«P2( 1)-P1( 1)>!'P2(3) >00003 800 
l + S(3)*{Pim*P2(2)-Pl (2)*P2( 1 ) ) 00003 810 


W = W1+.W2+W3 ' 00003820 

I F ( W • EQ«0 • ) GO TO 961 00003830 

W1 = Ul/W 00003840 

W2 = W2/W 00003850 

W3 = W3/W 00003860 

991 NCN(1,IN0D) = II 00003870 

NONIJ.INODl = 12 SLsO 

NGN(3.IN0C> = 13 00003890 

KCN(4,IN0C) = 14 00003900 



ooonr)Oooooon 


N0W(1,IN0D) = W1 
NCW(2, INOC) = W2 
NCWOtINOC) = W3 
NOW(A,INOC) = W4 

lF(UO.GT,0)WRITE(UOt202 )NODE{ INDDJ tNDDEK I 11 ,NOOET{ 12) ,NOOET( 13) , 
INCDEK 14) ,W1 tV^2,W3,W4 
ICOO CONTINUE ' 

RETURN 

END 


00003910 

00003920 

00003930 

00003940 

00003950 

C0003960 

00003970 

00003980 

C0003990 


SUBROUTINE TCOMP ( U I ».U01 , U02 1 IOT» NOD ».NODT , NON, NOW, NT I ME » T I ME , NODE ) 
COMPUTE AND OUTPUT STRUCTURAL NODE TEMPERATURES AT GIVEN TIMES. 
UI,U01,U02 = INPUT, OUTPUT FILE UNIT NUMBERS. 

IQT = FILE UNIT NUMBER FOR THERMAL NODE TJ ME-TEMPERA TLRE VECTORS. 
NOD = NUMBER OF STRUCTURAL NODES. 

NODT = NUMBER OF THERMAL NQDES- 

NGN(J,I) = JTH THERMAL NODE NUMBER (J=l-4) FOR STRUCTURAL NODE I. 
NOW(J,I) = JTH THERMAL NODE WEIGHT U=l-4) FOR STRUCTURAL NODE I. 
NTIME = NUMBER OF TIME VALUES FOR THERMAL NODES; 

TIME! I) = ITH TIME VALUE FOR THERMAL NODES. 

NCDE(I) = EXTERNAL I.D. FOR STR'UCTURAL NODE I. 

STOP 1001 = STRUCTURAL TIME IS OUTSIDE RANGE OF THERMAL TIMES, OR 
TIMES ARE NOT IN INCREASING ORDER, 

COMMON/COMTO/TO/COMTl/Tl 
C TCm,TllI) = TEMPERATURES OF ITH THERMAL NODE AT TIMES 0,1. 

INTEGER UI,U0l,U02, 1 0 T , NOD , NODT , NON ( 4 , I ) , NTI ME , NODE < 1 ) 

REAL N0W(4,1 ) ,TIME( 1 ) ,T0 ( 1 ),T1<.1) 

INTEGER I ,J,K,L,II ,ISTOP,IT,ISTOR(4) 

REAL STIME,C,W,TIM0,TIM1,F0,F1,ST0R(4) 

101 FCRMATIFIO.O, IlO) 

2Cl FORMATIF10i4,70(lH*) ) 

202 F0RMAT(4( I10,F10.:2) ) 

203 FORMATdH ) 

205 FORMAT(lHi) 

IF{U01.GT.O)WRITE(U01,205) 

REWIND lOT 
IT = 2 

READ! lOT) {Til D, 1 = 1, NODT) 

DO 2 1=1, NODT 


C0004000 

00004010 

00004020 

CC004030 

00004040 

00004050 

00004060 

00004070 

00004080 

00004090 

000C4100 

00004110 

00004120 

C0004130 

00004140 

00004150 

CC004160 

00004170 

00004180 

00004190 

00004200 

00004210 

00004220 

00004230 

00004240 

00004250 

00004260 

00004270 

00004280 
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2 TO ( n =Ti ( n 

READ( lOT) (Tl( n,I = l,NOOT) 

1 REAC(UI,101 >STIME» ISTQP 
IF( ISTOP.NE.O)GO TO 501 
IF{U01.GT«0)WRITE(U01,20l)STIME 
IF(U02.GT.0)WRITE(U02,201)STIME 
IF{STIME.LT.TIPE{IT-1 ) )STOP 1001 

3 IF { STIME.LE.TIMEl IT) )G0 TO 9 
IT = IT + 1 

IF( LT.GT-KTIKE)STOP 1001 
DO 6 I-=l»NODT 

6 Tom = Tim 

REAOnOT) (TUI ), I = l,.NQOT) 

GO TO 3 

9 T IMO = TIME! I T-I) 

TIKI = TI^EI IT } 

FO = TIMl-TIMO 

IFI FO.GT.Oi )F0 = (TIH1-STIME)/F0 
FI = l.O-FO 
I = 0 
11 L - 0 

DO 20 K=l,4 
I = I + l 

I F { I .GTiNQO )G0 TO 20 
L = L + 1 

ISTOR( L) = NODE ( I ) 

C = 0. 

DO 15 J=l,4 
I 1 =’ NONU, I ) 

W = NOW! J,I ) 

15 C = C + W*( F0*T0m ) + Fl*Tl( I I ),) 

STOR(L) = C 
20 CONTINUE , 

IF( L .NElO )G0 TO 51 

IF (UCl .GT«0 ) WRITE (UO 1,-203 ) 

IF(U02.GT,0) WRITE(U02,203) 

GO TO I 

51 IFCUOl.GT.O ) WRITE (UO 1,202 ) { I STOR ( K ) , STOR { K ) , K= 1 ,-L ) 
IF(U02.GT.0)WRITEIU02,202 ) ( I STOR( K ) „STOR ( K ) , K= 1 , L ) 
GO TO 11 
501 RETURN 
END 


00004290 

00004300 

00004310 

00004320 

00004330 

00004340 

0C004350 

00004360 

00004370 

00004380 

00004390 

00004400 

00004410 

00004420 

00004430 

00004440 

00004450 

00004460 

00004470 

00004480 

00004490 

00004500 

00004510 

00004520 

00004530 

00004540 

00004550 

00004560 

C0004570 

00004580 

00004590 

00004600 

00004610 

00004620 

C0004630 

00004640 

00004650 

00004660 

00004670 

00004680 

G0004690 

00004700 



